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We study the dynamics of two types of pairs of excitations which are bound despite their strong 
repulsive interaction. One corresponds to doubly occupied sites in one-dimensional Bose-Hubbard 
systems, the so-called doublons. The other is pairs of neighboring excited spins in anisotropic 
Heisenberg spin- 1/2 chains. We investigate the possibility of decay of the bound pairs due to 
resonant scattering by a defect or due to collisions of the pairs. We find that the amplitudes of the 
corresponding transitions are very small. This is a result of destructive quantum interference and 
Q I explains the stability of the bound pairs. 
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I. INTRODUCTION 



Cold atoms trapped in optical lattices offer high degree of control over such factors as the interaction strength, 
level of disorder, and system geometry. This makes it possible to simulate, in a controlled fashion, various models of 
condensed matter physics Bose and Fermi (fHH Hubbard models are among the studied ones, and spin-1/2 



systems are receiving increasing attention [9Hll|. Another significant advantage of cold atom systems is their weak 
coupling to the environment. This makes it possible to study their nearly dissipation-free evolution for long times, 
which in turn has led to the renewed interest in the far-from-equilibrium dynamics of many-body quantum systems. 

Several works have been recently dedicated to the formation [3! [i^ . [l3j |. dynamics (l^ - tl8| . and relaxation of 
bound pairs (BPs) of atoms that can exist even in the presence of repulsive interactions. They correspond to doubly 
occupied lattice sites and have been called doublons. Doublons have been observed in Bose Q and Fermi [l9[ 
Hubbard models with short range interaction. When the interaction is strong, be it attractive or repulsive, an energy 
gap can emerge between the band of doublons and the band where each atom is on a singly-occupied lattice sites. 
This gap, together with the isolation of the system, results in the long lifetime of the BPs in ideal systems. In some 
sense, this lifetime is even unexpectedly long. Similarly, BPs of parallel spins can emerge in spin-1/2 systems with 
strongly anisotropic nearest-neighbor interaction, commonly referred to as XXZ models. In a spin chain where all 
spins point down due to a strong magnetic field, for example, two up-spins on neighboring sites form a long-lived BP. 
£NJ \ All realistic systems have defects. They can lead to scattering of BPs. Understanding the stability of BPs in the 
presence of such scattering is therefore of central interest. For BPs in spin-1/2 chains with a single on-site defect and 
large anisotropic exchange interaction, this problem was addressed in our previous articles [2ll l22j. In this work we 
extend the analysis to doublons in one-dimensional Bose-Hubbard (BH) system and provide a comparison with the 
previous results. We also discuss the stability of BPs against many-particle scattering in ideal atomic and spin chains. 
The presence of an atom on a lattice site in the BH model or of a flipped spin in the XXZ chain is generally referred 
to as an excitation. 

We consider one-dimensional systems where the energy of the on-site positive Hubbard interaction or the nearest- 
neighbor S z — S z spin interaction largely exceeds the bandwidth of single-particle excitations. In the BH model, the 
bandwidth is determined by the hopping integral, whereas in the XXZ model it is determined by the in-plane spin 
coupling. This condition guarantees the energy separation of the BP band from the band of unbound two-particle 
excitations. We then analyze the case where the system has a site defect with an excess energy that also largely exceeds 
the bandwidth of single-particle excitations. Such a defect would be expected to cause strong scattering of BPs if a 
& • BP could resonantly decay into an excitation localized on the defect and a propagating excitation. This resonance 
can be achieved by tuning the defect site energy to a value close to the interaction energy. We show, however, that 
such inelastic scattering is largely suppressed. This is a consequence of destructive quantum interference. 

We also discuss the decay of BPs in a defect-free chain. The goal is to understand what keeps the BPs stable when 
they collide with each other Q • In the anisotropic Heisenberg model, energy conservation allows two colliding bound 
spin pairs to decay into a bound three-spin state and an unbound spin excitation. We show, however, that such 
scattering is also strongly suppressed by destructive quantum interference. For the BH model, energy conservation 
requires a multi-particle collision for the emergence of a bound triple of excitations. 
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This paper is organized as follows. The two models analyzed are described in Sec.|TTJ We consider strong interaction 
and a small number of excitations. The analytical studies developed here are based on Bethe ansatz [23-25] and 
perturbation theory, as explained in Sees. Hill and HV] The core of the work appears in Sees. |V] and EU In Sec. |V] 
we show that a BP has a parametrically small amplitude to split by sending an excitation to the defect site, and in 
Sec. I VII we show that a two-BP state essentially does not couple with a bound three-excitation state, even though 
both scattering processes are allowed by energy conservation. The BPs are almost elastically reflected from a defect 
site and from each other. 



II. MODELS 

We study one-dimensional systems with L sites and periodic boundary conditions (closed chains) described by the 
BH and XXZ models. In both cases we assume strong interaction and take H = l. 

A. Bose-Hubbard model 

In the BH model the interaction occurs on the same site and the Hamiltonian is given by, 

# BH = E ^ + if E n i ("j - !) - J E + b U h ) • 

Here, bj and bj are respectively bosonic annihilation and creation operators on site j, rij — b-bj is the corresponding 
number operator, J is the nearest-neighbor hopping integral, and U is the energy of the on-site particle interaction. The 
parameter e,- i s a site-dependent energy accounting for inhomogeneous external potentials, which can be experimentally 
controlled |26|. We assume that this background potential is the same for all sites, except for j = d, referred to as the 
defect site, so that Ej = e + g5j t d, with e > 0. The BH Hamiltonian conserves the total number of atoms (excitations) 
N, the dimension of each subspace being (L — 1 + N)l/[(L — 1)!AH]. 

B. XXZ model 

The anisotropic spin- 1/2 Heisenberg model, also known as the XXZ model, is described by the Hamiltonian, 

h xxz = E BjS] + j z E s j s Ui - J xy E i s i s u + S PU) • ( 2 ) 

3=1 j=l 3=1 

Here, Sj' v ' z = <jj' y ' z /2 are the spin operators on site j, a x .' v ' z being the Pauli matrices. The parameter Bj = B+gSj : d, 
with B > 0, gives the Zeeman splitting of each site, as determined by a static magnetic field in the z direction; it 
is the same for all sites, except for the defect site d. The Hamiltonian contains only nearest-neighbor exchange with 
coupling constants J z and J xy ; S*Sj +1 is the Ising-type interaction and SJSj +1 + is the flip-flop term [27| . 

We assume a large effective magnetic field, B> \J Z \ pointing up in the z direction, so that the ground state has all 
spins pointing down, independently of the sign of J z . A spin pointing up corresponds then to an excitation. 

The Heisenberg Hamiltonian conserves the total number of spins in the z direction, S z = J3t=i Sj< an d therefore it 
consists of blocks of dimension L1/\(L — N)\N\], where N is the number of up-spins. The defect-free system is intcgrablc 
and can be solved with the Bethe ansatz [23l - (25j [. It can also be mapped onto a system of spinless fermions [28| or 
hardcore bosons [2{|. In the presence of a single defect with excess ener gy q ~ J xy , J z , the system becomes chaotic, 
as attested by the Wigner- Dyson form of the level spacing distribution [30l - l32l ]. However, if the defect is very large, 
the system behaves as a chain with open boundary conditions, where the excitations move freely bouncing back and 
forth from the defect site. In this situation, analytical studies are again viable. 

In the next section, by combining perturbation theory and Bethe ansatz, we derive analytically the eigenvalues and 
eigenstates for the trivial subspace of a single excitation, where the BH and XXZ models are identical. It serves as 
an introduction for the sections where two excitations are considered. The basis vectors used are states where each 
excitation is on a single site. They are eigenstates of the interaction part of Hamiltonians ([TJ and (|2"j). as for example, 
|0010) for H ([]]) or equivalently | Uti) f° r H ©■ The table below gives their energies £n for the cases where the 
excitations are far from each other and away from the defect site. 
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TABLE I: Energies of the basis vectors for the BH and XXZ models, which correspond to -/V localized excitations far from each 
other and away from the defect site. 

III. ONE EXCITATION AND ONE DEFECT 

We start with a single excitation in a chain with a large-energy defect, \g\ >> \J\,\J xy \. The defect breaks 
translational symmetry. The scenario becomes similar to that of a system with open boundaries, the excitation either 
moves freely along the chain, being reflected when it gets next to d, or it is localized on the defect site. The results 
shown below are independent of the signs of g and J, J xy , so we assume them positive [2l[ |22j. 

A. A single freely propagating excitation 

For an excitation away from the defect, we look for approximate solutions of the type, 

d-l+L 

^F,free= E < X )<K X )> a(d) = a(d + L) = 0, (3) 

x=d+l 

where </>(x) indicates the basis vector where the excitation is on site x and a(x) is the probability amplitude of this 
state. 

The effect of the BH Hamiltonian (fT]) on 4>(x) is 

H Bii cf)(x) = e? u <j>(x) - J[<j>(x - 1) + cj>{x + 1)], (4) 

where £f H = s (cf. Table H|. Substituting Eqs. (|3]) and (j4]) into the Schrodinger equation Hipi = Eiipi and selecting 
the terms proportional to 4>(x), we obtain 

£f H a(x) - J[a(x - 1) + a(x + 1)] = E^a(x). (5) 

We now make an ansatz for a(x), taking into account the reflection of the excitation by the defect site, so that 

a(x) = Ce l6x + C'e~ lf>x , 

with the condition a(d) = a(d + L) = 0. This leads to the energies, 

E?? iee (6)=E? n -2Jcos0, (6) 

and the eigenstates, 

d-l + L 

^geeW= E A S m[6(x-dmx), (7) 

where A is a normalization factor and 

6 = — with k = 1,2,... (L- 1). 
L 

The same eigenstates are obtained for the XXZ model, the eigenvalues being 

O)^™- Jxy COSO, (8) 

with £^ xz given in Table |H 
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B. An excitation localized on the defect 



Since g is very large, the localization length of the excitation on the defect site is very small, which leads to the 
trivial state, V?£ c = "H^)' with energy 

^5c = fi BH +S+— • (9) 
9 

The last term above is a second-order perturbation theory correction, which when taken into account leads to excellent 
agreement with the numerics. It is derived from a virtual transition of the excitation on the defect to the site next to 
it and back to the defect site. For the XXZ model the energy of the localized state is given by 

J 2 

pXXZ _ cXXZ , n , u xy /, n x 

A i,ioc - fc i + 9 + t^-- (!UJ 
2# 



IV. TWO EXCITATIONS AND ONE DEFECT 



In the case of two excitations, we look for solutions of the type 



^2 BH -E a ( x ^)^y)' tit XZ = 52a(x,y)4>(x,v), (11) 

x<y %<y 

where <j)(x, y) indicates the state where the excitations are on site x and y and a(x, y) is its probability amplitude. In 
the BH model, there are states with x = y where two excitations are on the same site, while in the XXZ model, the 
excitations can reach neighboring sites, but cannot be on the same site, so x ^ y. 

We are interested in the case where the interaction is strong, \U\ J for the BH model and | J z \ 3> J xy for the XXZ 
model, which creates two well separated bands of energy. In one band, the excitations are noninteracting and move 
freely along the chain; in the other, the excitations are bound together due to the interaction and move slower as a 
single heavier particle [3, [U [22j ■ The cause of the onset of such pair of bound excitations is the difference in energy 
between the two bands. States with freely propagating excitations in the BH- model have energy ~ £^ H , whereas 
the BP-states have energy ~ £f K + U . For the XXZ model, the energies of these states are respectively ~ ^2 CXZ 
and ~£^ xz + J^. BPs can therefore emerge for either attractive or repulsive interactions. From now on, we assume 
repulsive interaction, U,J z >0. 

In addition to strong interaction, we consider also a defect with a large site energy, so that g 3> J, J X y and U — g 3> J 
for the BH model or equivalently J z — g ^S> J X y for the XXZ model. This leads to a third energy band placed between 
the previous two. Here, one excitation is trapped on the defect site and the other moves freely along the chain. The 
energy of this band is ~ £ 2 + <?. Thus, the three bands are separated in energy, as illustrated in Fig. [TJ The dashed 
lines in the diagram indicate two kinds of localized states that emerge close to the defect. These states and the energy 
bands are described in the next subsections. 



A. Two freely propagating excitations 

The lowest energy band in Fig.[T]is composed of states with noninteracting excitations; the excitations are sufficiently 
far apart from each other and from the defect site to move freely. In this case, the Schrodinger equation for the BH 
Hamiltonian leads to 

£i H a(x, y) - J[a(x -l,y) + a{x, y - 1) + a(x + 1, y) + a(x, y + 1)] = Eff iee a(x, y), (12) 

where £f H = 2e (cf. Table U). The defect works as a barrier that reflects the excitations before they reach site d. Due 
to energy conservation, the excitations are also reflected before reaching the same site. An appropriate ansatz for the 
probability amplitude away from the defect is then 

a(x, y) = Ce i9x e- i9 ' y + C e~ if,x e ie ' y , (13) 

with the condition a(d, y) = a{x, d + L) = a(x, x) = 0. By substituting (|13|) back into (|12j) . we obtain the eigenvalues 

*%free(0> 6 ') = £ 2 BH - 2J(COS0 + COS0'), (14) 
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FIG. 1: Two- excitation energy spectrum for the BH model when U ~S> J and U/2 < g < U. We assume g,J,U> 0. The band 
centered at £™ = 2e is formed by noninteracting excitations. The band centered at 2e + g + 2J 2 /g is formed by states with 
one excitation localized on the defect site and the other moving freely. The narrow band centered at 2e + U + 4 J 2 /U is formed 
by propagating pairs of bound excitations. The dashed line above the BP-band indicates the localized states where the BPs 
are next to the defect; E BUI f = 2e + U + 2J 2 U/[g(U — g)] (see text). The dashed line at highest energy indicates the BP on the 
defect site. 



constituting a band of width 8 J. 

There are (L — 1)(L — 2)/2 states with two freely propagating excitations for the BH model and (L — 2)(L — 3)/2 
states for the Heisenberg model. For the latter, the energies are 

£*free(M') = £ 2 XXZ - J xy (cos0 + cosO'). (15) 



B. An excitation localized on the defect and a freely propagating excitation 



The energy of the states of the BH model where one excitation is localized on the defect and the other moves 
freely between sites d + 1 and d + L — 1 (referred to here as loc-free-states) differ by g from the states where both 
excitations propagate freely and by U — g from BPs. By taking into account the energy of the excitation on the defect 
up to second order of perturbation theory, the eigenvalues and eigenvectors of the loc-free-band can be found from 
the equation 



£ 2 BH + .9 



J 2 f„ r U 2 -g 2 + 4Ug \ 



i(d, y) - J[a(d, y-l) + a{d, y + 1)] = E^ oc _ frcc a{d, y) 



To the lowest order in J 2 /g, J 2 U/(U 2 — g 2 ), we seek the solution in the form 

a(d,y) = Ce l6y + C'e- l9y , 



(16) 



(17) 



with the condition a(d,d) — a(d,d + L) = 0. The scenario is very similar to that of a free excitation described in 
Sec. IIIII The eigenvalues are well approximated by 



^2,loc-free(0) 



£ 9 BH 



2J 2 



2Jcos6> 



and the eigenstates have the form 



?/; BH 
r2,loc-frcc 



L 



d-l+L 

(6)= Asm[0(y~dMd,y), 
y=d+i 

with k = 1,2,... (L- 1), 



(18) 

(19) 
(20) 
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where A is a normalization factor. 

The results for the Heisenberg model are similar [see Table [TT| for a mapping between the two models], except that 
here there are only L — 3 states, since the boundary conditions are a(d, d + 1) = a(d, d — 1 + L) = 0. 

C. Extended and localized bound pairs 

In the leading order of perturbation theory, transport of a pair of bound excitations through the chain occurs via 
an intermediate virtual transition where the pair splits and then recombines again. The BP behaves as a single heavy 
free excitation in an open chain. In the BH model, it moves from site d+1 to d + L — 1 with effective hopping integral 
2J 2 /U. The eigenvalues and eigenvectors can be found from the equations 



rBH 
C 2 



rBH 
& 2 " 



.-BH 



-U 

■U -\ 
■U -\ 



4J 2 
f — 
4J 2 

— 4 
4J 2 

— 4 



2J 2 

IT 



a(x, x) H — ^-[a(x — l,x — 1) + a(x + 1, x + 1)] = .Ef B P a(x, x), 



2J 2 



(21) 
(22) 



i{d + 1, d + 1) + ^a(d + 2, d + 2) = £^i P a(d + 1, d + 1), 
2 7 2 

(d+L-l,d+£-l)+— a(d+L-2,d+L-2) = Ef$ P a(d+L-l,d+L-l), (23) 



where 



2J 2 



2J 2 



2J 2 ff 
C/((7- ff ) 



is the energy difference between the states right next to the defect, 4>{d + 1, d + 1) and </>(d + L — 1, d + L — 1), and 
the other BPs. This difference arises because, in the second order of perturbation theory, the virtual dissociation of 
a BP with one excitation hopping onto a defect site gives a different energy denominator than in the case where the 
virtual transition is made onto a regular site. 
The ansatz 



a(x, x) = Ce l9x + C'e 



in Eq. (|21"T) leads to the eigenvalues 



d < x < d + L, 



4J 2 

IT 



4J 2 

IT 



(24) 



They form a narrow energy band of width 8 J 2 /U. 

Equation (|2~lj) has to be compatible with Eq. (|22|) for x = d + 1 and with Eq. ([23]) for x = L + d — 1 . This leads to a 
system of two linear equations for the coefficients C and C . The equations have a nontrivial solution for 9 satisfying 
the condition 



f{6) - /(-*), f(0) = e^ L ~V 



2J 2 

IT 



l 2 



(25) 



As the defect energy increases, two states eventually split off the band of extended BP-states. They correspond 
to surface states, where the excitations are localized primarily next to the defect, on sites d+1 and d + L — 1. In 
the thermodynamic limit, L — > oo, the energies of these states, E SUT {, are derived from the two imaginary roots of 
Eq. (USD, 



Im 6 = ±ln 



fJ 



U-g 



which appear for U > g > U/2. These energies are shown in Fig. [T] with a dashed line. We then have two surface 
states and L — 3 extended BPs where the pairs move between sites d + 2 and d — 2 + L. A similar scenario emerges 
for the XXZ model. When J z > g > J z /2 there are two states with excitations localized next to the defect, on sites 
(d + 1, d + 2) and (d + L — 2, d + L — 1), and L — 4 extended BP eigenstates. 

A localized state of higher energy appears when a BP is placed on the defect site. For the BH model, this is a single 
state with energy £f K + U + 2g. For the Heisenberg system there are two such states with energy £* xz + J z + g; 



7 



they have excitations on sites (d,d + 1) and on (d, d + L — 1). These states can form symmetric and antisymmetric 
linear superpositions [2l| HU . 

Table \H\ summarizes the results for both models for the loc- free-band and the band of extended BPs. These states 
together with the BPs localized next to the defect are the main parties to the quantum interference effects analyzed 
in Sec. El 



BH model XXZ model 



U > J and U/2 < g <U J z » J xy and J z /2 < g 


< Jz 


One excitation localized on the defect and one freely propagating excitation 


£ 2 B L_ frco = £ 2 BH + g + 2J 2 /g- 2J cos 9 E™ z c _ ilcc = £ xxz + g + J 2 y /(2g) - J xy 


cos 9 


Extended bound-pairs of excitations 


£ 2 B Ip = £ 2 BH + U + 4J 2 /U + (4J 2 /U) cos Ef^l = £ 2 XXZ + J z + J 2 xy /{2J Z ) + \J 2 xy /{2J z )] 


cos 9 



TABLE II: Energy spectrum of loc-free and extended BP-states in a chain with a defect and two excitations. The energies of 
basis vectors with two separated excitations away from the defect for the Bose-Hubbard and the anisotropic Heisenberg models, 
£ BH and £ xxz , respectively, are given in Tabled] 



V. ANTIRESONANCE STABILIZATION OF BOUND PAIRS IN A CHAIN WITH A DEFECT 



In physical systems of interest for the studies of doublons, the defect energy can often be tuned. These systems 
display a somewhat counterintuitive behavior if the energy of a loc-free state resonates with the energy of a BP away 
from the defect, that is if g is close to U in the BH-model. In this case the band of extended BPs lies inside the (much 
broader) band of loc-free states. One could then expect hybridization of the two types of states. In other words, a BP 
could approach the defect and resonantly scatter into an excitation localized on the defect and a propagating single 
excitation, or vice versa, a propagating single excitation could resonantly scatter off the excitation on the defect, and 
then the two excitations would move away from the defect as a bound pair. Both scenarios are allowed by energy 
conservation, but they have very small probabilities to occur, because of a destructive quantum interference. 

When g w U ^> J, BPs localized immediately next to the defect are directly coupled to the loc-free-states. They 
form a new hybrid band denoted below with the subscript "hyb" . The eigenstates and eigenvalues are obtained from 
the equations, 



[£ 2 BH + g] a(d, y) - J[a(d, y - 1) + a(d, y + 1)] = E™ yh a(d, y), 
[£jp H + U] a(d + 1, d + 1) - V2Ja(d, d+ 1) = E^ yh a(d - 
[£jp H + U] a(d+L-l,d+L-l)-V2Ja(d+L-l,d) 



with a{d, d) = a(d + L,d + L) = 0. By manipulating the ansatz 



^2, hyb 



l,d + l), 
,a(d+L-l,d+L-l), 



a(x,y) 



Ce i6{y+1) + c / e -ie(y+l) 



(26) 
(27) 
(28) 

(29) 



\ r< id(d+l) i ril — X X U i6(d+L+l) , W -iOld+L + l)] c r 



with d+l<y<d+L — 1, in the equations above, we derive the relation between the coefficients C and C, 



c 



2i8d 



C 



(U - g)e lB + J 
{U - g)e- lf> + J 



(30) 



and the energies 



£2 B hyb=£ 2 BH +S-2JC0S< 



(31) 



with 9 obtained from 



m = /(- 



f(9) = e WL \U-g + Je 



,01 
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In the case of exact resonance, g = U, the eigenstates are approximately those given by Eqs. (|19|) and ([20]) . but now 
with the sum in y from d+ltod + L + 1 and k = 1, 2, . . . L + 1. 

We might expect the BP-states to mix with the states from the hybrid band above in second order of perturbation 
theory through the states <j)(d + 2, d + 2) and <p(d + L — 2,d + L — 2). The BP on site d + 2 [and similarly for site 
d + L — 2] would either hop onto site d + 1 or split into sites (d, d + 2). However, the probability amplitudes for 
these transitions are almost equal in magnitude and opposite in sign, therefore canceling out. This can be seen from 
Eqs. (HU) and ([30]), which lead to 

2iC sin(29)e iSd 

a(d + l,d + l) + a(d,d + 2) = — — [U - g + 2Jcos0] . (32) 

{U g)e + J 

The sum above approaches zero when U « 5 — 2 J cos to zeroth order in as used to obtain Eq. (|3"Tj) . The absence 
of hybridization is therefore a result of destructive quantum interference and guarantees the stability of the BPs. A 
similar anti- resonance occurs also for BPs in the XXZ model [U, [22| . 

One can think of the phenomenon as a Fano anti-resonance, where the narrow band of BPs in the presence of a 
defect goes through the broad band of loc-free states without scattering. This absence of scattering can be equivalently 
understood from the following analysis of the motion of the bound pair on site d + 2 (the case of the site d + L — 2 is 
completely similar). The Schrodinger equation for this pair is 



£ BH + jj + ^1 _ E BH 

u 



2 J 2 

a(d + 2, d + 2) - —a{d + 3, d + 3) - s/2Ja(d + 1, d + 2) = 0, (33) 



where E BB w £ BB + U. Above, we have kept the term of hopping to site d + 3, as in Eq. (|2T|) . but not the hopping 
to d+ 1; instead, we have the amplitude for the intermediate virtual state <j){d + 1, d + 2). Even though this last term 
is nonresonant, it is directly coupled to the following resonant terms, 

(£f H - E BB ) a(d + 1, d + 2) - V2J [a(d + 1, d + 1) + a(d + 2, d + 2)] - Ja{d, d + 2) = 0. (34) 

Notice that a(d + 1, d + 3) is nonresonant. 



We now substitute Eq. (|34]l into Eq. ([33]), taking into account that E BH — £f H « U. We also use the Schrodinger 
equations for the states near the defect, with energies E BIi , 

[£f H + g - £ BH ] a(d, d + 1) - \/2 Ja(d + 1, d + 1) - Ja(d, d+ 2)] = 0, (35) 
[£ 2 BH + £/ - £ BH ] a(d + 1, tZ + 1) - V2 Ja(d, d + 1) = 0, (36) 

to rewrite a(d, d + 2) in terms of a(d + 1, d + 1). [Eq.(27) and (36) coincide if one sets £; BH = EW b ]. This finally 



leads to 



U 



2 J 2 

a(d + 2, d + 2) - —a(d + 3, d + 3) 



w -(J- 1 [£ 2 BH + C/ - £ BH ] [£ 2 BH + .g - £ BH ] o(d + 1, d + 1). (37) 

For propagating BPs, the energy E BH lies within the band centered at £ BH + U + AJ 2 /U with width 8J 2 /U (see 
Fig. rjj, therefore the coefficient at a(d + l,d+ 1) in Eq. ([37]) is - (J 2 /U)\g - U\/U < J 2 /C7. This means that the 
propagating bound pairs do not hop onto site d + 1 (nor onto site d + L — 1). They are elastically reflected from the 
defect and do not hybridize with the loc-free states. 

Figure [2] illustrates the decoupling of the resonating BP- and loc- free-bands for the BH model (top panels) and the 
XXZ model (bottom panels). A small chain is considered to facilitate the visualization of the states involved. The 
figure shows the time evolution of the probability P to find a state <j)(x, y) in the full two-excitation wave function 
of the system ^(t), with P = P(x,y;t) — | (<j)(x, y)\^f(t)) \ 2 . The initial states on the left panels are BPs away from 
the defect site, whereas for the middle and right panels they are a loc-free-state. On the left, only BP states are 
seen, confirming the stability of the pair of bound excitations even when its splitting is not prohibited by energy 
conservation. The pairs move back and forth through the chain, being reflected by the defect. The middle and right 
panels show that the states that can display an appreciable coupling with the initial state are of two kinds: (i) the 
states with one excitation on the defect site and the other moving back and forth in the chain (middle panels), and 
(ii) the BPs localized next to the defect (right panels). 



9 



Initial state: BP 



Initial state: loc-free 





1 

0.8 



\(1,3) A( 15 ) 






i_(i,4)/ \ 


- (2,3) 


(5,6) 




"/. V. \. 





100 200 
J t 



300 



2 



AT 



4 6 
J t 

xy 



8 10 



FIG. 2: Time evolution of two-particle excitations in the BH model (top) and XXZ model (bottom) for a chain of length L — 6 
with a defect on site d = 1. The curves give the probability P to find a state <f>(x,y) at time t. The defect excess energy is 
exactly equal to the interaction strength: g — U for the Bose-Hubbard model and g — J z for the XXZ model. The repulsive 
interaction is strong, U/J = J z /Jx V = 30, as in the experiments of Ref. |4j. Left panels: the initial states are bound pairs 
(doublons) with wave functions 0(3,3) (top) and 0(3,4) (bottom). Middle and right panels: the initial state is a loc-free state 
with wave function 0(1, 3) (top and bottom). The bound-pair states on the right panel are localized on sites next to the defect 
("surface states"), no BPs in the bulk have an appreciable amplitude. 



VI. ANTIRESONANCE STABILIZATION OF BOUND PAIRS IN A DEFECT-FREE SYSTEM WITH 

MANY EXCITATIONS 

We now consider more than one BP in a defect-free chain, g = 0. The goal is to analyze whether the pairs remain 
stable when colliding or can form larger clusters of bound excitations. We show that, similarly to the previous 
section, destructive quantum interference prevents BPs from merging into bound triples of excitations, even when 
this transition is energetically favorable, and therefore, as long as their density is not too high, BPs remain long-lived 
many-particle states. 

Of particular interest for not too high densities is the decay of BPs as a result of a collision of two BPs. In the 
BH-model with strong interaction, U ^> J, the energy of a state with two BPs is ps 2U. Since the total number of 
excitations is conserved, the two BPs could decay into a bound triple of excitations and a freely propagating excitation 
or in a bound quartet of excitations. The energy of a triple of excitations on the same site is « 3U and the energy 
of a quartet is ps 6U. Therefore, by energy conservation, a collision of two BPs (two doublons) cannot lead to their 
decay in this model. Decay processes would require a collision of at least 3 doublons. In the XXZ model, on the other 
hand, the decay of BPs could indeed result from a collision of two BPs. In this case, the energies of the bound triple 
(3 upward spins on neighboring sites) and of two BPs are ps 2J z . In principle, the two BPs could then combine into 
a bound triple and a free excitation. The amplitde for this transition is, however, very small, as we show below. 

To study the stability of the BPs in an ideal XXZ chain, we present a four-excitation wave function as 

1pf XZ — ^ a(xi,X2,X3,X4)(f>(xi,X2,X3,X4) ) (38) 

where </>(xi , x%, X3, x±) is the state of the system with spins on sites x±, ...,xn pointing up, whereas all other spins are 
pointing down. An eigenstate with two spatially separated bound pairs is described by the wave function (|38p where 
the coefficients have the structure a(x, x + 1, y, y + 1) with y > x + 3. These states form a narrow band centered at 
rj £f xz + 2J Z with typical width - 2J% y /J z , cf. Table II. 

A bound triple and an unbound free excitation are described by the wave function (|38[) with coefficients that have 
the structure a(x, x + 1, x + 2,y) with y > x + 4. They also form a band centered at ps £ xxz + 2J Z , but with a much 
larger width ~ J xyi which is determined primarily by the hopping of the free excitation. 

Two BPs separated by only one site are directly coupled to a triple. To the leading order in J xy / J z , the Schrodinger 
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equation for the corresponding resonating amplitudes for energies E reads 

(£f xz + 2J Z - E xxz ) a(x, x + 1, x + 2, x + 4) - -y^[a(x, x + 1, x + 3, x + 4) + a(x, x + 1, x + 2, x + 5)] = 0. (39) 

Here, the triple is on sites (x, x + 1, x + 2) and the BPs are on sites (x, x + 1) and (x + 3, x + 4). 

BPs separated by more than one site could mix with triples and with BPs separated by one site in second order of 
perturbation theory, but the amplitudes of these processes turn out to be very small. To describe what happens, we 
consider the BPs on sites (x, x + 1, x + 4, x + 5) and show that they have a very small amplitude of hopping onto sites 
(x, x + 1, x + 3, x + 4) [or equivalently onto site (x + 1, x + 2, x + 4, x + 5)], that is, into the state strongly admixed 
with bound triples. To this end we write the Schrodinger equation 

( J 2 \ J 2 

£f xz + 2J z + ^f- E xxz a(x, x + 1, x + 4, x + 5) + -ff- [a(x, x + 1, x + 5, x + 6) + a(x - 1, x, x + 4, x + 5)] 

-- J X j/ [a(x, a; + 1, a; + 3, x + 5) + a(z, x + 2,x + i,x + 5)] =0. (40) 

The last two terms in this equation are nonresonant and correspond to a virtual dissociation of one of the BPs. They 
describe the intermediate states that lead to hopping of the BPs toward each other. The Schrodinger equation for 
a(x, x + 1, x + 3, x + 5), taking into account Eq. (|3T))) and the fact that E xxz « £ xxz + 2J Z for two-BP states, leads to 

^J xy a(x, x + 1, x + 3, x + 5) w (f^ xz + 2 J z - £ xxz ) a(aj, a; + 1, x + 2, x + 4). (41) 

Function a(x, x + 2, x + 4, x + 5) is similarly expressed in terms of the bound-triple amplitude a(x + l,a;-l-3,x + 4, a; + 5). 

It is immediately seen from Eqs. (|40|) and (|41l) that, for the characteristic energies E xxz lying within the bandwidth 
of the BP excitations, \E XXZ — (£ xxz + 2J Z )\ < J 2 y /J z , the admixture of the states with BPs separated by two sites 
and the states with bound triples is parametrically smaller than the bandwidth of two-BP-states. Physically, this 
means that once two BPs approach each other, they are reflected before they become closer than two empty lattice 
sites. Such reflection is due to a destructive quantum interference between the bound triple and two-BP states 
separated by one site and can also be interpreted in terms of Fano anti-resonance. 




FIG. 3: Time evolution of 4-excitation states in a defect-free XXZ chain of length L — 9 for large repulsive interaction 
Jz/Jxy =30. The curves show the probabilities P — (4>(xi, X2, £3, X4) |V'? XZ (' : )) | 2 f° r different sites. The initial state 0(1, 2, 6, 7) 
couples only with states from the two-BP-band, so only states with BPs separated by two or more sites are seen. 

The time evolution of the four-excitation states in the XXZ model with a large repulsive interaction is illustrated 
in Figure [3] for a defect-free chain with nine sites and periodic boundary conditions. In this system, there are nine 
states with two BPs separated by two or more sites. The initial state considered has the BPs separated by three sites, 
^(1,2,6,7). These BPs move along the chain and reflect before the inter-pair distance becomes equal to 1. In the 
figure, only the nine curves corresponding to the states of the two-BP-band participate in the evolution, confirming 
the decoupling between this band and the hybrid band of triples and BPs separated by one site. Not a single state 
with three excitations next to each other or with two BPs separated by one site had an appreciable amplitude in the 
simulation for the studied time range. This is in full agreement with the analytical results. 



VII. CONCLUSION 



We studied the dynamics of bound pairs of excitations in one-dimensional systems with strong repulsive short-range 
interaction described by the Bose-Hubbard and the anisotropic Heisenberg models. In the first system, the bound pair 
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is formed by two atoms occupying the same site. Such state has become known as doublon. In the second system the 
bound pair corresponds to parallel spins placed on neighboring sites. For strong interaction, the bandwidth of bound 
pairs in an ideal chain is much smaller than the bandwidth of states with free excitations, so the bands of unbound 
and bound excitations are well separated. 

Of primary interest was the question of stability of the bound pairs in two possible decay scenarios. One refers to a 
chain with one bound pair and a single on-site defect. The energy of the defect can be tuned close to the interaction 
energy of the bound pair so as to allow energy conservation in the scattering of the BP off the defect. The second 
scenario refers to an ideal chain where the total energy of two BPs is equal to the sum of the energy of a bound triple 
of excitations and the energy of a propagating single excitation, so that we might expect the BPs to merge into a 
triple. 

Our central result is the observation that, for strong interaction, both processes are suppressed. This is because 
there are two intermediate states that can mediate the decay process of the pairs and these states have essentially 
equal and opposite in sign amplitudes. Such destructive quantum interference leads to the long-time stability of the 
pairs of bound excitations. 

The stability of bound pairs of excitations with strong negative binding energy has immediate bearing on the 
possibility of observing them in the experiment. In experiments with cold atoms it is possible both to create a 
"defect" in the optical lattice and to control its energy. It is also possible to change the density of bound pairs. An 
observation of the stability of the pairs would be a direct demonstration of quantum antiresonance in the system. 
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